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' Thermonuclear supernovae are valuable for cosmology but their physics is not yet 

fully understood. Modeling the development and propagation of nuclear flame is 
complicated by numerous instabilities. The predictions of supernova light curves 
^ \ still involve some simplifying assumptions, but one can use the comparison of 

the computed fluxes with observations to constrain the explosion mechanism. In 
spite of great progress in recent years, a number of issues remains unsolved both 



■ in flame physics and light curve modeling. 
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p ,■ 1 Introduction 

O 

^ \ Supernovae of type la (SNe la) are important for cosmology (better to say, 

^ \ for cosmography) due to their brightness. They are not standard candles, but 

can be used for measuring distances with the help of the peak luminosity - 
decline rate correlation, established by Yu.P. Pskovskii |^ and M.M. Phillips 



pOl (see the review 0). To exclude systematic effects in linking the observed 
5^ \ light of distant SNe la to the parameters of cosmological models, one has to 

understand the nature of supernova outbursts and to build accurate algorithms 
for predicting their emission. 

This involves: identifying the progenitors of SNe la; the birth of thermonu- 
clear flame and its accelerated propagation leading to explosion; light curve 
and spectra modeling. 

When the full understanding will be achieved, one can try to evaluate the 
importance of evolution effects in using supernovae as distance indicators. In 
spite of great progress in recent years, a number of issues remains unsolved 
both in flame physics and light curve modeling. We point out some problems 
which seem most important to us. 
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2 Progenitors 



Mechanical equilibrium and evolution of stars is easily understood from the 
virial theorem for a star: 3 J PdV = —U , where P is the pressure, V is volume, 
and U is the gravitational energy of a star. Crude estimates V ~ i?^, M ~ 

R^p, U GM'^/R in the virial equation give P - GM'^/^p^/'^ . For an 

ideal classical gas, the equation of state P = TZpT implies TZT ~ GM'^^^p^^^, 
with TZ the gas constant. This is already enough to understand the evolution of 
massive stars! While a star loses energy and contracts, its internal temperature 
T grows. If the losses are balanced by a nuclear energy release, then the 
contraction stops and a thermal equilibrium is established: 

nuclear heating power L"^ = radiative cooling (luminosity) L~. 

The rate of thermonuclear heating scales as {avo) ~ exp[— (ac/T)^/^] due to 
the Gamow's peak: the chances to penetrate the Coulomb barrier for fast 
nuclei grow, but the tail of Maxwell distribution goes down. Here ac depends 
strongly on nuclei charges Zf. ac oc thus high-^ ions can fuse only 

at high T. Small perturbations of T produce huge variations in L"*" since, 
normally, T ^ ac- 

The cooling L~ depends on T moderately, and it seems that L~ cannot 
compensate an overheating perturbation. Why then do not all stars explode 
violently? The reason is the same as for the growth of T in stars losing the 
energy at contraction: they have negative heat capacity. For a star made of a 
classical plasma with 7 = 5/3 the internal energy Q = \nMT and the virial 
theorem implies U = —2Q, so the total energy £ is negative: £ = U + Q — 
—Q < 0. Thus any growth in £ due to heating leads to the drop of T (because 
nuclear energy is used for expansion against gravity) . The perturbations decay, 
and there is no hope to get a thermonuclear supernova from a normal star 
composed of a classical plasma. 

The situation changes, if a star is made of a degenerate matter. Then the 
heat ^TZMT resides in ions, but its absolute value is much less than Q which 
is now governed by Fermi energy of electrons. We have crudely for a mixture 
of non-relativistic electron Fermi gas and classical ions: P ~ Kp^/^ + TZpT ~ 
GM'^/^p^/^, and TZT ~ GM'^'^ p^/^ - K p^'^ . So at high density, the equilibrium 
temperature T decreases with growing p, i.e. goes the same way as £. The 
total heat capacity becomes positive, and runaway can set in as in terrestrial 
explosives. So, a progenitor of SN la must be a degenerate star - a white 
dwarf. 

A single white dwarf is unable to explode, it cools down. But when it is in a 
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binary system the chances to produce a supernova do appear (we need only one 
in ~ 300 dying white dwarfs to explode in order to explain the rate of SNe la). 
Even if the binary has two dead white dwarfs, it can explode because they can 
merge due to emission of gravitation waves (double-degenerate, or DD scenario 
|1TI| ). If one star in the binary is alive, the white dwarf can accrete its lost 
mass and reach an instability (single-degenerate, or SD scenario |27, |^). It is 
unclear which scenario is most important, there are strong arguments |p.4|| from 
chemical evolution that only SD is the viable one. On the other hand, it seems 
that DD can produce a richer variety of SN la events. Moreover, discoveries 
of intergalactic SNe la 0, can be explained more naturally, because a DD 
system may evaporate from a galaxy. It is quite likely that both scenarios are 
being played, but their relative role may change in young and old galaxies. 
If so, a systematic trend may appear in SNe la properties with the age of 
Universe, and this may have important consequences for cosmology. 



3 Thermonuclear flames 

After merging in DD scenario, or after the white dwarf accretes large amount 
of material in SD case, the explosive instability develops. In principle, com- 
bustion can propagate either in the form of a supersonic detonation [|I| wave. 



or as a subsonic deflagration [|T2|, (flame). In detonation, the unburned fuel 
is ignited by a shock front propagating ahead of the burning zone itself. In 
deflagration, the ignition is governed by heat and active reactant transport, 
i.e. by thermal conduction and diffusion. 

Most likely, the runaway starts as a laminar flame propagating due to ther- 
mal conduction. In terrestrial flames, the 'fusion' of molecules goes with the 
rate: (avo) ~ exp{—Ea/TZT), - the Arrhenius law of chemical burning. Here 
Ea is activation energy. The parameter, showing the strong T-dependence of 
the heating Ze = (91og(o"fo)/i91ogT ~ Ea/TZT is called the Zeldovich num- 
ber in the theory of chemical flames. For them typically Ze ~ 10 . . . 20. The 
classical theory |^ predicts the flame speed V{ ~ Ze~^[vTlT /Trea.c{Tb)Y^'^ , with 
Tre3.c{T) oc exp[Ea/ {'R-T)]. In SNe, for nuclear flames, rreac(^) oc exp[a^'''^/ (3T-'^/^)], 
and, Ze = 51og(crfo)/91ogT ~ a^^ /{3T^^^), which has values very similar to 
terrestrial chemical flames. 

A big difference with chemical flames is the ratio of heat conduction and 
mass diffusion, the Lewis number, Le = {vtIt) / {vdId)- One finds Le ~ 1 in 
laboratory gaseous flames, while Le ~ 10^ in thermonuclear SNe, since heat is 
transported by relativistic electrons, vt ~ c, and there is almost no diffusion. 
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Table 1: Flame speed V{ and width Zf in C+0 



p 




h 


Ap/p 


Ma 


10^ gcc 


km/s 


cm 






6 


214 


1.8 X 10-5 


0.10 


2 X 10-2 


1 


36 


2.9 X 10-* 


0.19 


4 X 10-3 


0.1 


2.3 


2.7 X 10-2 


0.43 


4 X 10-* 



It S> /z). Nevertheless, the modern computations follow the old theory ^0 



closely. The conductive flame propagates in a presupernova with ff which is 
too slow to produce an energetic explosion: the ratio of V{ to sound speed, i.e. 
the Mach number. Ma, is very small (see Table 1). The star has enough time 
to expand, to cool down, and the burning dies completely. So an acceleration 
of the flame is necessary in order to explain the SN phenomenon. This is the 
main problem in current research of SNe la hydrodynamics. 

There is a rich variety of instabilities that can severely distort the shape of 
a laminar flame. The Rayleigh-Taylor (RT) instability governs the corrugation 
of the front on the largest scales. On the smallest scales the flame is controlled 
by the Landau-Darrieus (LD) instability. RT, LD instabilities and turbulence 
make computations difficult, but without them a star would not explode. All 
these instabilities were considered already by L. Landau [|l^ as a means to 
accelerate the flame. 

Because of instabilities, the flame surface becomes wrinkled and its area 
grows as S* oc -R" , with average radius R and a > 2, i.e. faster than S oc R^. 
In other words the surface becomes 'fractal'. The exponent a is actually 
the fractal dimension, a = Dp. The effective flame speed is determined ^8 



by the ratio of the maximum scale of the instability to the minimum one: 
Ves = (Amax/Amin)'^''"^- When the vorticity is not important it is possible to 
study in detail the non-linear stage of LD instability and to find the fractal 
dimension [Q. A similar dependence of the flame fractal dimension on the 
density jump across the front was found in SPH simulations of the flame 
subject to RT instability |^. 

The fractal description is good for LD while it remains mild, because it 
operates in a star on the scales from the flame thickness (a tiny fraction of a 
cm) up to ~ 1 km. For the RT instability, Amax/Amin is very uncertain and 
the fractal dimension is uncertain too. So a direct 3D numerical simulation is 
necessary. The same is true for a low density regime of LD when it is strongly 
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coupled to turbulence (generated on the front, or cascading from large RT 
vortices). A great progress is achieved here in several groups 



0, m, y]- 

When simulating 3D turbulent deflagrations one encounters two problems: the 
representation of the thin moving surface separating hot and cold material, 
and the prescription of the local velocity ff of this surface as a function of the 
large-scale flow with a crude numerical resolution > 1 km. One solution is 
sketched in 



22 1; for a different approach see [13|. In spite of the progress this 



problem cannot be treated as completely solved, and even ID approach may 
give interesting results, especially for unusual SNe la [^. 



4 Light curves of SNe la 

Given a hydrodynamic structure of SN ejecta, one can compute a light curve 
which should be compared with observations. There are several effects in SNe 
physics which lead to difficulties in the light curve modeling of any type of SNe. 
For instance, an account should be taken correctly for deposition of gamma 
photons produced in decays of radioactive isotopes, mostly ^^Ni and ^^Co. To 
find this one has to solve the transfer equation for gamma photons together 
with hydrodynamical equations. Full system of equations should involve also 
radiative transfer equations in the range from soft X-rays to infrared for the 
expanding medium. There are millions of spectral lines that form SN spectra, 
and it is not a trivial problem to find a convenient way how to treat them even 
in the static case. The expansion makes the problem much more difficult to 
solve: hundreds or even thousands of lines give their input into emission and 
absorption at each frequency. 

Currently, powerful codes appear aimed to attack a full 3D time-dependent 
problem of SN la light Yet there are some basic questions, like averaging 
the line opacity in expanding media, that remain controversial. 

In our work we predict the broad-band UBVI and bolometric light curves 
of SNe la, using our ID-hydro code which models multi-group time-dependent 
non-equilibrium radiative transfer inside SN ejecta. In our previous analysis 
we have studied two Chandrasekhar-mass models: the classical defiagration 
model W7 ]18|] and the delayed detonation one DD4 |2^, as well as two sub- 



Chandrasekhar-mass models: helium detonation model LA4 11161 and low-mass 



detonation model WD065 with low ^^Ni production p^, which was constructed 



for modelling subluminous SNe la events, such as SN 1991bg. All those models 
were simplified spherically-symmetrical (ID) ones. 

The UBVI light curves of ID models are shown in Fig. |T|. The Chandrasekhar- 
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Figure 1: UBVI light curves for 4 ID models (W7 - solid line, DD4 -dots, LA4 - 
dashes, WD065 - dash-dots). Crosses, stars and triangles show the light curves for 
three observed SNe la. 



mass models demonstrate almost identical light curves. The sub-Chandrasekhar- 
mass ones are more different. WD065 has almost similar element distribution 
as Chandrasekhar-mass models, and the shape of its light curve is in principle 
the same as that of W7 and DD4. It is just much dimmer due to a very low 
^^Ni abundance. 

LA4 is very different from any other model, since the explosion there started 
on the surface of a white dwarf, not in the center, as for every other model, so 
there is a ^^Ni layer near the surface in LA4. This feature explains why the 
model is essentially bluer than other ones. 



Currently we employ our new corrected treatment for line opacity [24 



the expanding medium, which is important especially in UV and IR bands. 
It seems that ID thermonuclear supernova models, e.g. the deflagration W7 
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Table 2: Parameters of SN la models 



Model 


DD4 


W7 


LA4 


WD065 


MR 


Mwd" 


1.3861 


1.3775 


0.8678 


0.6500 


1.4 




0.63 


0.60 


0.47 


0.05 


0.42 




1.23 


1.20 


1.15 


0.56 


0.46 



in Mf 



^'in 10^^ ergs s ^ 



18 model and the delayed detonation DD4 [29] one, produce the light curves 



fitting the observations not so good as the recent 3D deflagration model MR 
computed at MPA P2[. We believe that the main feature of the latter model 



which allows us to get the correct flux during the first month, is strong mixing 
that moves the material enriched with radioactive ^^Ni to the outermost layers 
of SN ejecta. 

Fig. ^ demonstrates that in spite of quite different structure of the old W7 
model and the new MR one their light curves are similar in many details. 
Moreover the new model behaves better in U and B bands. Unfortunately, 
the bolometric light curve for MR model is somewhat too slow. The ejecta 
must expand with a higher speed to let photons to diffuse out faster. 



5 Conclusions 

There are several points which require attention for applying SNe la in cos- 
mology: progenitors may be different in younger galaxies; burning regimes 



may change with the age of Universe [^. The physical understanding of the 
Pskovskii-Phillips is not yet achieved. 



In the flame modeling the new 3D SN la model ||22[ is very appealing. Yet 



it is not a flnal one: a detailed post-processing of nucleosynthesis is not yet 
checked in the light curve calculation. 

The SN light curve modeling still has a lot of physics to be added, such as a 
3D time-dependent radiative transfer, including as much as possible of NLTE 



effects ||T0|, which are especially essential for SNe la. All this will improve our 
understanding of thermonuclear supernovae and their role in cosmology. 

Our work is partly funded by RFBR (grants 00-02-17230 and 02-02-16500). 
SB is grateful to Ana Mourao for support. 
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Figure 2: UBVI light curves for the 3D (MR; sohd) and ID (W7; dashed) models. 
Crosses, stars and triangles show the light curves for three observed SNe la. 
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